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ABSTRACT 

We discuss a new method which is potentially capable of constraining cosmological 
parameters using observations of giant luminous arcs in rich X-ray clusters of galaxies. 
The mass profile and the mass normalization of the lenses are determined from X-ray 
measurements. The method also allows to probe the amount and equation of state of 
the dark energy in the universe. The analysis of a preliminary sample of 6 luminous, 
relatively relaxed clusters of galaxies strongly favours an accelerating expansion of the 
universe. Under the assumption that the dark energy is in the form of a cosmological 
constant, the data provide an estimate of f^Ao = 1-1 with a statistical error of ±0.2. 
Including the constraint of a flat universe and an equation of state for the dark energy 
wx ^ —1, we obtain fluo — O.lOiO.lO and wx — — 0.84±0.14. Relaxing the prior on 
Wx, we find that the null energy condition is satisfied at the 3-a confidence level. 
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, A growing body of evidence suggests that most of the energy 
' density in the universe consists of some sort of dark energy. 
' Recent observations of the power spectrum of fluctuations 
in the cosmic microwave backg round (CM B; ISper gel et all 
i2003h 1. of large sca le structure iHaw kins et al as well 

as supernovae data jKnop et all200^]^nr^^^Ll20o3l de- 
mand for significant pressure-negative dark energy that is 
not clustered with the galaxies. 

The phenomenological equation of state wx = px/px, 
i.e. the ratio between the pressure px of the un- 
known energy component and its rest energy density 
px has therefore assumed a central role in observa- 
tional cosmology. The simplest explanation for dark en- 
ergy is a cosmological constant, for which wx = —1; 
another possibility is quintessence, i.e., a scalar field 
rolling down an almost fiat potential fCaldwell et al. '199^; 
[Ratra and Peebles 1998; Rubano and Scudcllaro 2002); k- 
essenc e, namely a scalar field with non canonical kinetic 
terms iArmendariz-Picon et al]ll999^ . or models based on 
branes and extra dimensions, such as the Cardassian sce- 
nario llZhu and Fuiimotdliool) . can also drive an acceler- 
ated expansion. 

Even though the equation of state is generally time- 
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dependent, a properly calculated weighted mean of wx can 
ap proximate the late time behavior of a large class of mod- 
els iWang et al ."2000^). In a Friedmann-Lemaitre- Robertson- 
Walker (FLRW) model of universe, the distance depends 
on Wx only thro ugh a multiple integral on the redshift 
jMaor et alJ200ll) . and current data do not enable to recon- 
struct the time evolution of the equation of state. In what 
follows, we will consider an equation of state for the dark 
energy assumed to be constant with redshift. 

Investigations on wx are usually restricted to a limited 
range by assuming that the so called null energy condition 
(NEC) of general relativity is fulfilled on macroscopic scales. 
Under quite general hypotheses, NEC pu ts a lower bound 
to the value of lOx llSchuecker et alj|2003l) . While no strict 
mathematical proofs of NEC exist, from a phenomenological 
point of view this assumption is not justified. 

The necessity of extending dark energy analyses is 
further illustrated by some models of quintessence with 
a time dependent wx{t) ^ —1 in which the assumption 
of constant wx leads to an estim ated value of the mea n 
equation of state smaller than — 1 iMelchiorri et al.ll2003f) . 
So, it is useful to enlarge the parameter space to wx < 
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Matter with wx < —1 is named phantom energy 
JCaldwell et alj|2003l) . Theoretical motivations for such an 
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unusual behavior have aheady been propos ed, cf. for in- 
stance the string theory dCaldwell et al.ir2003l) . The positive 
phantom energy density becomes infinite in a finite time, 
overcoming all other forms of matter and leading to a big 
rip. 

Although several independent methods of observations 
all converge towards a coherent picture, the real signature 
of the accelerated expansion of the universe is provided by 
the apparent magnitude versus redshift test for type la su- 
pernovae. Other attempts to build the Hubble diagram, 
like those based on the compact radio source angular size 
(jChen and Ratra 2003), give constraints that, even though 
consistent with supernovae data, are much weaker. System- 
atic effects, such as luminosity evolution, gray intergalactic 
dust, gravitational lensing or selection biases, could mimic 
the effects of dark energy llSereno et alJ 120021: Ebnrv et alJ 
L2OO3). Although it is believed that these contaminants are 
under control, an independent constraint can improve the 
statistical significance of the statement about the expansion 
of the universe: it is therefore still useful to develop new 
tools for the determination of the cosmological parameters. 
Gravitational lensing studies can provide such methods. 

Results from the statistical strong gravitational lensing 
of flat-spect rum radio sources based on the CLASS sample 
JChadl2003l) are in good agreement with those from type 
la supernovae. Here, we are interested in cluster of galaxies 
acting as lenses. Clusters of galaxies provide a laboratory 
for studying and measuring the energy content of the uni- 
verse in a variety of ways. They are the largest virialized 
structures and their importance in observational cosmology 
is well known. In the beginning of the last century, the need 
for unseen dark matte r was first stat ed by direct estimates 
of their total masses llZwickvlll933ll . After a century, as- 
tronomers face again darkness. The turn of dark energy has 
came. 

Clusters of galaxies acting as lenses on background 
high redshift galaxies have been lon g proposed as sources 



nign redsnitt galaxies nave been lon g proposed as sources 
of in fo rmation on the uni ve rse jPac z^ak^mi Gorskil 
19811: iBreimer a nd Sanders' '1992^ 'Fort et al.l | l997|: 
Onk^an ^ 1 9 98; Lombardi an d B erth] ^93 : 
Gautret et al.l I2OO0I : iMeneehetti et airBoojh . Values of 



a particular ratio of angular diameter distances can be 
determined once the modeling of the lens is constrained 
and both th e arc p o sition and i t s reds hif t are m e asured 
JChiba et al.l l2002l: iGolse et alj |20o2). ISerend (120021) 
discussed the feasibility of using these gravitational lens 
systems to discriminate between accelerating or decelerating 
models of universe and to probe the amount and equation 
of state of the dark energy. The observations of a suitable 
number of lensing clusters at intermediate redshift can 
determine the equation of state. 

The combined analysis of independent observations of 
clusters of galaxies at different wavelengths can greatly 
improve the knowledge of the system allowing to gain 
an insight into the sec o nd-order cosmo logical parameters 
JPe Filippis et alJ 120041) . ISerenol J2003D used a combined 
analysis of the intracluster medium (ICM) within the hydro- 
static region, as derived from both Sunyaev-Zel'dovich effect 
measurements and X-ray images, to constrain the Hubble 
constant and the cosmological matter density. 

Under suitable hypotheses, in a cluster of galaxy, X- 
ray and lensing observations probe the same potential well. 



Here, we propose the use of X-ray spectroscopic and surface 
brightness analyses of clusters of galaxies to determine the 
mass profile and the mass normalization. In Sec.|21 we show 
how the mass profile of a cluster can be derived from X-ray 
measurements. Section |3 discusses some threshold require- 
ments for strong lensing events in X-ray clusters. In Sec. 01 
we present the cluster data sample used in the statistical 
analysis and the selection criteria applied. Section |3 is de- 
voted to deriving cosmological parameters and testing the 
NEC condition. A discussion on some systematics that can 
affect the method is contained in Sec.|^ Section [7| contains 
some final considerations. 



2 ISOTHERMAL /3-MODEL 

At least in theory and under some simple assumptions, the 
distribution of the cluster total mass can be inferred from 
the modeled gas pressure distribution. In fact, the ICM is 
usually assumed to be isothermal and in hydrostatic equi- 
librium in this potential, whereas nonthermal processes do 
not contribute significantly to the gas pressure. Under the 
assumption of spherical symmetry, we have 



GM{r) _ ksTx d 
fiTUp dr 



In nc(r-). 



(1) 



where M(r) is the total mass of a cluster within radius r; 
Tx and rieir) are the gas temperature and number density, 
respectively; rrip is the proton mass; /x denotes the mean 
molecular weight, assumed to be constant throughout the 
gas, so that the electron number density traces the gas den- 
sity. Its value is usually fixed with solar metallicity mea- 
surements and we shall therefore assume ^ = 0.585, with an 

error of 4%. 

A conv entional /3-model dCavaliere and Fusco-Femiand 
Il976lll978h . which is widely adopted to fit the density profile 
of ICM, is: 



no(r) = rioo 1 + 



-3/3x/2 



(2) 



where noo is the central electron number density, Tc the core 
radius and the parameter /3x determines the slope. By as- 
suming a constant X-ray spectral emissivity, Eq. ((21 predicts 
the following X-ray surface brightness, 



6x(6l) =bxo 1 + 



-3fe-l-l/2 



(3) 



which provides a good fit to X-ray observations; 6 (= r/D^), 
where Dd is the angular diameter distance to the cluster, is 
a dimensionless angular variable. 

Inserting Eq. (j^J in Eq. (0, we get the total mass of 
the isothermal /3-model, 



M{r) = 



3fcBTx 



The mass density distribution is given by, 

W = 2 = P07 ^^^^> 

An r'' dr f, / \2l 



1 + 



(4) 



(5) 



where 
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pa = 



3fcB Tx/^x 



(6) 



4TvGfj,mp r'i 

The surface mass density, projected along the hne of sight, 
is 

V 2 



3/2 ' 



where Eq is the central surface density, 
^ 3 fca Tx/9x 1 

^0 = X 



2 Gfinip dc -Dd 
The mass enclosed by a cylinder of radius 8 is, 



Mcyi(6») = 27rEo (ScDd) 



1+ f 



1/2 ' 



(7) 



(8) 



(9) 



The profile in Eq. ^ is known as softened power-law 
surface mass density a nd the slope parameter is fixed to 1/2 
^Schneider et aUFmol. 

In a sphericall y symmetric, regular le ns, a strong lensing 
event can occur if dSchneider et al.(ll992ll : 



Eo > Ecr, 

where Ecr is the critical surface mass density. 



(10) 



(11) 



being Dds the angular diameter distance from the lens to 
the source and Da the angular diameter distance from the 
obser ver to the source. A t angential critical curve appears 

at et ^^^^EmI^^, 



1/2 



a radial critical curve forms at 



(ecJ 



2/3 



1 1/2 



(12) 



(13) 



The corresponding radial caustic in the source plane has an 
angular radius rjr, where 



(14) 



A source inside the radial caustic produces three images, one 
which instead is outside has only one image. 

The mass profile of a cluster can be fully determined 
by X-ray measurements. A morphological analysis of the X- 
ray surface brightness allows to constrain the parameters Oc 
and /3x, while the gas temperature Tx is determined with 
spectroscopic measurements. 



3 STRONG LENSING CRITERIA 

As seen in the previous section, when we assume a spherical 
mass distribution, a lens can produce multiple image sys- 
tems only if the central surface density overcomes the criti- 
cal one. According to Eqs. l(8l lll|l . the criterium in Eq. I|10|l 
is fulfilled if 
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Figure 1. The ratio of distances Ds/D^^ as a function of wx in a 
fiat model of Universe with Qmo = 0-2. Here, we assume = 0.3 
and 2= = 1. 
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Figure 2. The ratio of distances Ds/Dds as a function of Qmo in 
a fiat model of Universe with cosmological constant (uix = — 1). 
Here, we assume = 0.3 and Zg = 1. 
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In other words, as it was to be expected, very hot and mas- 
sive clusters with a small core radius are more likely to form 
arc-like systems. 

A proper X-ray modelling is needed to check the con- 
dition in Eq. (1151 . Estimating the temperature profile re- 
quires a detailed spatially resolved spectroscopy. A multi- 
phase structure in the central ICM could result in significant 
discrepancies between emission-weighted and mass- weighted 
estimates of the temperature. A simple isothermal model 
which does not incorporate the cooling-flow components or 
other complexities in the X-ray structure can lead to sys- 
tematic underestimates of order 10-40% in the temperature 
and mass determinations (Allen 1998). 

On the average, detailed X-ray morphologies show rapid 
variations in the central regions; such variations may be due, 
for instance, to recent merger or accretion activities and, 
therefore, they might invalidate the assumption of hydro- 
static equilibrium. Complex dynamical activities can lead 
to overestimating the core radius of the dominant clump in 
the cluster, as inferred from X-ray data. The core radius 
for the gravitating matter in cooling-flow clusters is usually 
small, in agreement with optical studies. The radius nearly 
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doubles for intermediate systems, and, in the case of non - 
cooling- flows, it may be larger by a factor six iAllerJll99d) . 

Threshold values for strong lensing events depend on 
the cosmological parameters through the ratio of distances 
Ds/Dds- We consider a FLRW universe filled with non- 
interacting pr essureless matter and dark energy, describe d 
as X-matter JChiba et al.lll997l: iTurner and Whit3ll997l) . 
Since the contribution from relativistic particles is negligi- 
ble in the redshift range investigated in our analysis, we will 
neglect it in what follows. In such a model of universe, the 
angular diameter distance between an observer at Zd and a 
source at Zs is 



D{zd, Zs) 

where 

£{z) = 



1 



1 



rSinn 



£{z) 



■dz 



(16) 



H{z) 
Ho 



(17) 



= Vf^M0(l + Z)3 + nxo(l + 2)3(»x + l) + J7ko(1 + 2)2, 

and Ho is the present value of the Hubble parameter; f^MO 
and Qxo are the today normalized densities of dust and dark 
energy, respectively; Qko = 1 — flmo — Hxo; Sinn is defined 
as being sinh when Qko > 0, sin when Qko < 0, and as the 
identity when Qko = 0. For the expressi on of the distance in 
an in homogeneous universe, we refer to lSereno et alJ l)20nil 
|2003). 

Depending on the values of the cosmological parame- 
ters, significant variations in the distance ratio have to be ex- 
pected. As it can be seen from Figs. the ratio Ds/Dds 
is an increasing function of both Qmo and wx- Dark energy 
with large negative pressure helps the formation of arc-like 
systems in X-ray clusters. The best case is that of a cosmo- 
logical constant. Changing wx from —1 to 0, the ratio of 
distances increases of 13% for a lens at Zd = 0.3 and a back- 
ground source at = 1 in a fiat universe with Qmo =0.2. 
Sub-critical values of the pressureless matter density also 
favour strong lensing systems. The relative variation from 
f^MO = to r^MO = 1 is ~ 21% for wx ~ —1. A large value of 
a cosmological constant increases the cross section for strong 
lensing and lowers the threshold value of the combination of 
X-ray parameters in Eq. HSU . 



4 CLUSTER SAMPLE 

Our method to determine cosmological parameters relies on 
a number of assumptions: hydrostatic equilibrium, isother- 
mality, spherical symmetry. The selection of the sample 
must therefore be very careful in order not to violate them 
|SmitheEai]|2003). We have collected data from literature 
on luminous, relatively relaxed clusters with giant luminous 
arcs and well measured. X-ray spatially-resolved observa- 
tions. Then, we have selected a sub-sample suitable for our 
analysis, by requiring: i) good agreement of optical and/or 
lensing mass measurements with the X-ray results; ii) a 
regular X-ray morphology; in) temperature measured with 
Chandra Observatory or XMM-Newton satellites; iv) arcs 
with known spectroscopic redshift. 

The agreement between mass estimates from indepen- 
dent methods firmly limits the systematic uncertainties 
when determining the cosmological parameters. Apparently 



relaxed clusters can sho w significant discrepancy between X - 
ray and lensing masses (iKneib et alJl2003l: lOta et al liooi). 

Together with the regularity of the X-ray surface brightness 
emission, the consistency among the X-ray, optical and/or 
gravitational (both weak and strong) lensing mass measure- 
ments is a reliable indication of hydrostatic equilibrium and 
strongly suggests that the global properties of the clusters 
can be used as cosmological probes. X-ray gas should be 
supported by thermal pressure which dominates over non- 
thermal processes and the assumption of hydrostatic equi- 
librium can be considered valid. 

Gas temperatures are typically very hard to mea- 
sure (iGioia et alj|2004) . whereas low resolution measure- 
ments suffer from point-source emission contaminating clus- 
ter spectra, Chandra's or XMM-Newtori's spatial resolutions 
allow the identification of point sources and an independent 
spect ral de termination of the temperature in each radial 
sheU llAUen et al.ll2001al: iDonahue et alJl2003l) . In order to 
be conservative in our selection procedure, we consider only 
clusters for which the temperatures determined with last 
generation satellites agree with previously reported values 
from ROSAT and/or ASCA data. Such an agreement is a 
strong indication that the temperature estimate is substan- 
tially free from systematics and that the adopted uncertainty 
is realistic. 

We wish to stress that conditions i) and ii) are often 
fulfilled by cooling flow clusters, which turn out to be quite 
regular, dynamically relaxed systems. The X-ray emission 
is typically symmetric and little or no substructures appear 
at optical wavelengths. These clusters are also amongst the 
most X-ray luminous. 

Furthermore, the comparison of masses derived from X- 
ray and stron g lensing show excelle nt agreement for cooling- 
flow clusters llAllenlll99dlwJ2000h . while discrepancies ap- 
pear in dynamically active clusters with either small or no 
cooling-flows. The oversimplification of a spherical model 
and the application of isothermality and equilibrium hy- 
potheses may be inappropriate in presence of local dynami- 
cal activities. 

Due to their regularity, from the point of view of gravi- 
tational lensing, cooling-fiow clusters are "critical", i.e. they 
may form giant arcs, only if the condition in Eq. (1151 is ful- 
filled. On the other hand, substructures and non-equilibrium 
states in dynamically active, non-cooling flow systems en- 
hance the probability of detecti ng strong lensing events 
jBartelmann and Steinmetd 1 19961) , which are favourite in 
the presence of irregular and complex dark matter distri- 
butions. A simple spherical model, which does not account 
for the local dynamical activities, may lead to overestimat- 
ing the gravitating masses inferred from lensing studies. 

The nearly coincidence of a X-ray single peak with the 
cD galax;y position, which has been argued to be a good 
indicator of the presence of local violent activities of the in- 
tracluster gas in the central core, has als o been prop osed as 
an indicator of hydrostatic equilibrium llAllenlll998t) . Typi- 
cally, in cooling-flow clusters the emission is sharply peaked 
to a position nearly coincident or very close to the optically 
dominant cluster galaxy, which generally is the centre of the 
matter distribution as inferred from lensing studies. 

After the application of our selection criteria, our fl- 
nal sample reduced to 6 clusters where strong lensing anal- 
ysis and X-ray studies are likely to trace the same grav- 
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Table 1. Information on arcs and X-ray prop erties of the strong lensing cluster sample. " References: (]j_[^lei^^aljj222^|) ; (2) 

t Schmidt ct al. 12004); (3) Donahue ct al. (2003) : (4) lArabadiis ct al, (,2002.1 : (5) Allen et al. (2001b.') : (&) iMolikawa et aI.lTl999l) : (7) 
Hicks et al]il2002l) : fS) lAUen et all l2002^ . 



Cluster 


M 


^arc 




T(KeV) 


/3 




ref" 


A 2390 


0.228 


0.913 


38.0 




0.48±0.02 


11.8±0.3 


1,2 


MS 0451.6-0305 


0.539 


2.911 


25.9 


10.3+I;g 


0.70±0.02 


31±1.0 


3 


MS 1358.4+6245 


0.328 


4.92 


21.0 


7.16±0.10 


0.69±0.02 


20.5±5.2 


4 


MS 2137.3-2353 
PKS 0745-191 
RX J1347.5-1145 


0.313 
0.103 
0.451 


1.501 
0.433 
0.806 


15.3 
18.2 
34.9 


5 5R+0'46 
M.OOq -,5 

12.2±0.6 


0.47±0.02 
0.57±0.02 


18.01°;^ 
4.8±0.3 


5,6 
5,7 
8 



itational potential well. The main properties of this sam- 
ple are listed in Table Q Our sample, spanning the range 
0.1 < Zd < 0.54j constitutes an updated sub-sample of the 
list in IWul (l2000l) . In formation on positions and redshifts of 
arcs were taken from lWul i200(]|) and references therein. New 
spectroscopic redshift measurements for MS 0451.6-0305 
and M S 2 137.3-2353 were in stead taken from iBorvs et alJ 
J2004 and lSand et~ai] J2002t) . respectively. 

For each cluster, the ICM has been modeled with an 
isothermal /3-model. We consider fits to the X-ray surface 
brightness which extend well beyond the core radius. When 
more morphological analyses are available for a single clus- 
ter, we adopt results from fits that extend to the outer ra- 
dius. Since we collected data from independent groups, in 
order to avoid unrealistic small uncertainties when compar- 
ing the data, we assume an error of 0.02 on the slope /3 when 
determined from Chandra observations. 

Details on the X-ray spectroscopic and imaging anal- 
ysis for each cluster can be found in the references listed 
m Table El Clusters in our sample have been the subjects 
of X-ray, optical and, often, Sunyaev-Zel'dovich effect stud- 
ies. They present very regular X-ray isophotal patterns. De- 
spite some local substructures which have been revealed in 
highly detailed X-ray surface brightness maps, morphologi- 
cal analyses support the spherical approximation: there are 
no evidences for violent merger shocks and the luminos- 
ity and tempe ratures do not seem to h a ve been boosted 
substantially [mi on ot al.' '2001a', '2002*: Arabadiis ct alJ 
b002; D onahue et aL200a : .Hicks ct al..2002:.Molikawa et al. 

Apart from MS 0451.6-0305, all the clusters in our final 
sample present massive cooling fiows in the central regions. 
Furthermore, for each cluster except again MS 0451.6-0305, 
X-ray morphological analyses have revealed a sharp cen- 
tral surface brightness peak at a position nearly coincident, 
within the astrometric resolution of Chandra observations, 
with the optical centroid for the dominant cluster galaxy. 
In MS 0451.6-0305, whereas the soft-energy peak of the X- 
ray surface brightness lies on t he brightest cluster ga laxy, a 
harder peak is slightly shifted dDonahue et al ]!2003). How- 
ever, due to the consistency between the X-ray and Sunyaev- 
Zel'dovich estimates of the central electron density and the 
agreement between the X-ray, len sing and optical estim ates 
of the total mass of this cluster iDonahue et aljl2003l) . we 
have included it in our final sample. 
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Figure 3. Contours of constant Ax^ in the Omo-^^AO plane. Lines 
are traced for Ax^ values of 2.30, 6.17 and 11.8. Accelerated mod- 
els of universe are above the dashed line. The full line represents 
the locus of flat models of universe (f2K = 0); the long-dashed 
line represents models with a null distance to Zs = 4.92; bouncing 
models in the upper-left shaded region do not have big bang. 

5 DATA ANALYSIS 

Let us now perform a statistics. Observations of giant 
arcs in X-ray clusters enable us to estimate the distance 
ratio DdJDs. From Eqs. l|Hl[niIIlJ, we get 

The parameters Tx, /3x and 9^ are derived from a X-ray data 
analysis of temperature and surface brightness profiles. 
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Figure 4. Confidence contours on the nMO-w'x plane. The shaded 
regions show the joint 1, 2 and 3-(t confidence regions, bounded 
by Ax^ values of 2.30 and 6.17 and 11.8, respectively. Accelerated 
models of universe are below the dashed line. 



The value of the angular radius of the tangential criti- 
cal curve was derived from t he observed arc radius Sarc • Ac- 
cording to lOno et al.l il999h . a conventional analysis where 
Ot = Obic, could overestimate the strong lensing mass by 
10-30%. To correct for this effect, we take 9t — eOaic, with 
e= ±0.04. 



Thex 




XO , WX ) 



(19) 



The statistical errors ai include uncertainties arising from 
Tx, /3x, dc, fJ, and Sarc- Many input data in our analysis are 
presented with asymmetric uncertainties. To obtain unbi- 
ased estimates, we apply correction formulae fo r mean and 
standa rd deviation as given by Eqs. (15, 16) in iD'Agostiiiil 
i2004f) . Since the analysis of the X-ray parameters is highly 
correlated, we add the uncertainties from /3x and 9c as max- 
imum errors. 

Once we assume that the description of the data 
is a good approximation, the likelihood of the cosmo- 
logical parameters turns out to be jC{Qmo , ^xo , wx) oc 
g-x^(nMo."xo.™x)/2 Following Bayesian statistics, the pos- 
terior probability of a particular model after analyzing the 
available data is proportional to the prior probability of that 
model multiplied by the likelihood. For each parameter in 
the model, its one dimensional posterior probability func- 
tion is found by marginalizing over all other parameters. 
Unless otherwise noted, for each parameter we quote the 
one-dimensional mean and standard deviation. 

We first consider the case of a cosmological constant 
{Qxo = f^Ao, wx = — 1). We limit our analysis to the region 
in the parameter space for which the distance to the most 
distant arc in our sample (at Zs = 4.92) is not null. In Fig.|21 



we show the function. According to the maximum likeli- 
hood ratio method, joint 1, 2 and 3-cr confidence regions for 
two parameters are bounded by Ax^ values of 2.30 and 6.17 
and 11.8, respectively. As it can be seen, our method is quite 
sensitive to the cosmological constant, while the pressure- 
less matter density is poorly constrained. The best-fitting 
model, (nMO,f2Ao) — (0.01,0.99), is nea rly coincident wit h 
a de Sitter model. As it was discussed in lGott et al.l i200ll) . 
appropriate priors should be as vague or noninformative as 
possible. We consider two types of priors: i) a uniform prior, 
with the prior probability proportional to di}ModQAo', H) 
the logarithmic or Jefi'rey prior, appropriate for a positive 
variable like JImo, with a prior probability proportional to 
df2AodfiMo/f^MO- Results on Qmo are highly sensitive on the 
choice of the prior, indicating that bette r data are neede d 
to convincingly constrain this parameter llGott et al.ll200lll . 
On the other hand, the determination of the cosmological 
constant is nearly insensitive to the prior. 

Let us consider an uniform prior. Within 3-a statistical 
significance, the data strongly favour accelerating models, 
whereas an Einstein-de Sitter model of universe {Qmo = 
1, JIko ~ 0) is not consistent. After marginalization, we find 
^Ao = 1.1±0.2. 

Then, following theoretical prejudices on an inflation- 
ary, nearly flat model of universe with a positive cosmo- 
logical constant, we restrict our analysis to the square re- 
gion of the parameter space defined by S5 f^MO ^ 1 
and ^5 Qao ^ 1. Adopting a uniform prior, we find 
Qmo = 0.2±0.2 and Oao = 0.87±0.12. Whichever prior we 
use, our analysis provides a strong evidence supporting the 
existence of dark energy. 

Let us now analyse the dark energy equation of state. 
We consider fiat models (SIk = 0). As a first step, we im- 
pose a prior on the equation of state, — 1 ^ wx ^ 0. The best 
fit model occur for {Q,mo,wx) — (0.01,-1). Adopting uni- 
form priors, after marginalization, we get Qmo = O.lOiO.lO 
and wx = — 0.84±0.14; accelerating models of universe are 
favoured at the 2-a level, see Fig. 01 Our analysis favours a 
model of universe in accelerate expansion with the energy 
budget dominated by dark energy. At this point, with the 
prior on wx, the cosmological constant turns out to be the 
best candidate as dark energy. 

Our results are in a substantial agreemen t with 
current determinations fHannestad and MortscU ' ' 



iMelchiorri ct al. 


2003; KnoR et al. 2003; SDcrgcl ct al. 


200| 


ISchuecker et al.l l2003l: iTonrv et al.l 


20031 


IZhu et al.1 l2004r) 



A combined analysis of supernovae, galaxy distortion 
from the Two-Degree Field Galaxy Redshift Survey (2dF- 
GRS; Hawkins et alj J2003M and CMB data from WMAP 



(jSpergcl et aL 2003 i) provide an estimate of Qmo = 0.27ln'n^ 
and a 95% upper confidence limit of wx < 0.78 (Knoc 



l2003h . An independent analysis in lSoergel et al.l 

some different external constraints and combining CMB, 
2dFGRS power spectrum, supernovae, and the HST Key 
Project data, also gives the same constraint on wx- 



5.1 Energy conditions violations 

Under quite general hypotheses, in the case of a spatially 
fiat FLRW geometry, with a negligible contribution from 
relativistic parti cles, NEC puts a lower bound on wx, wx ^ 
-(1 - Qmo)~'^ dSchuecker et alJl2003ll . For f^Mo = 0, we 
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Figure 5. Confidence contours on the f^MO-^x plane. Shaded 
regions are as in Fig. |1] The dot indicates the best fit model. The 
NEC condition is fulfilled above the dashed line. 

have wx 5S —1, i.e., the usual limit assumed in most of the 
investigations. 

Following ISchuecker 6^11 J2003l) . we use lux to test 
NEC. We find a best fitting value of tux ^ —3.28. Whereas 
at 2-G confidence level the NEC is violated, at the 3-cr con- 
fidence limit, it is fulfilled, see Fig. |^ The case of a cosmo- 
logical constant is consistent with data at the S-cr level. 

Distance data alone do not well bound from below the 
equation of state and the lower limit of tux is quite uncertain 
( Knop ct al. 2003). The weight of probability at very low 
u;x pulls the confidence level downward. From supernovae 
data, at 99% confidence level, an upper limit of wx < —0.64 
is de rived, but interval extends to wx < — 10 iKnop et alJ 
|20o3). The degeneracy can be broken by using orthogo- 
nal constraints, such as those arisi ng from abundances o f 
nearby clusters of X-ray c l usters dSchuecker et alj l2003l) . 
The analysis in lKnop et al.l ll2003fl provides a measurement 
of wx = —l.OSlg'jQ (statistical) ± 0.09 (identified syst emat- 
ics). The independent analysis in lSoergel et alJ i2003r) gives 
a similar result, wx = —0.98 ± 0.12. 



6 SYSTEMATICS 

In the previous section, we have performed a statistical anal- 
ysis based on data found in the literature. Whereas we se- 
lected a regular sample, with well measured properties, a 
number of systematics can still plague our analysis. We want 
now to address how they can affect our results. 

6.1 Density profile 

We have our high-redshift sample by using a /3-profile. Appli- 
cation of such a model provides robus t estimates of the total 
mass, as shown for simulated clusters jMohr et alll999^ and 



observed clusters jEttori et al.ll20o3) . It is still interesting, 
however, to investigate the impact of a central cusp in the 
dark matter distribution, as predicted by numerical simula- 
tions, on our estimates. To this aim, we consider the so called 
universal Navarro, F renk and White (NEW) density profile 
iNavarro et al.lll995fl . Since the similarity between the dis- 
tribution of ICM tracing the dark halo of the NEW poten- 

^^^^^^^^^^^^^^^^^^^^^^ 1f^^^^^^^ ^^^^^^^^^^^^^^"^^^^^^^^^^^ 

tial and the /3-profile iMakino et al. 1998:' Ettori fc EabianI 
[T9991I . NEW parameters can be derived from the /?- model 
ones with helpful emp irical formulae . We t herefore applied 
the relations found in iMakino et alj il99d) . The introduc- 
tion of a central cusp slightly reduces the estimated value of 
the cosmological constant of ~ 0.07 and the best fit value of 
Tiix increases of ~ 0.3. 

6.2 Temperature 

Isothermal models provide a good re presentation of th e 
physical properties of galaxy clusters iEttori et aljl200^ . 
Whereas the presence of a gradient in the temperature pro- 
file would reduce the total mass measurements and increase, 
consequently, the derived value of the cosmological constant, 
expected gradients in clusters are not st eep enough to affec t 
significantly the estimated total mass llEttori et al.ll2002ll . 
However, in order to be more confident on the temperature 
values used in our analysis, we checked our sample for consis- 
tency of the new temperature measurements obtained with 
high resolution satellite with those previously reported from 
ROSAT and/or ASCA data, covering a larger field of view 
(see data in Table 2 in Wu (2000)). This should control and 
exclude any systematic bias. In any case, a systematic offset 
in the temperature estimate of ~ 5% would determine a vari- 
ation in Q.AO of ~ 0.15, with a positive offset determining an 
overestimation of the cosmological constant. The constraint 
on lux, when no priors are adopted, would be nearly washed 
out by such an offset. 

6.3 Ellipticity and substructures 

One of the more important sources of indetermination comes 
from the modeling of the mass profile of the lens. Morpho- 
logical X-ray analyses revealed a tendency for clusters of 
galaxies to be elliptical, with a mean projected axis ra - 
tio of ~ 1.25 jMohr et alJ Il995l : iDe Eihppis et all 1200411 . 
Furthermore, together with the overall mass profile, sub- 
structures should also be considered. High detailed X-ray 
surface brightness maps can reveal interesting details in the 
distribution of ICM. In RX J1347.5-1145 |Allcn ot al. 2002), 
Chandra observations revealed evidences of shocked gas to 
the south-east of the main X-ray peak, coincident with a re- 
gion of enhanced Sunyaev-Zel'dovich effect, probably result- 
ing from recent subclustor merger activity. Presence of some 
subs tructure in A 239 (Allen ct al. 2001a) and PKS 0745- 
191 jHicks et alJl2002h was also detected and a hint that the 
brightest cluster galaxy may be just settling into the cluster 
core has been reported for MS 0451.6-0305 jPonahue et alJ 
l2003h . 

However, when the lens presents a rather regular mor- 
phology, even if a "not correct" potential shape is used in the 
reconstruction, or the contribution of small sub-structures 
is neglected, the cosmological paramet ers are still retriev ed, 
although with somewhat larger errors JColse et alJl2002h . 
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To best investigate the role played by these deviations 
from the assumed projected cylindrical symmetry, we con- 
sidered an additional source of systematic uncertainty in the 
estimation of the position of the tangential critical radius, 
of the order of ~ 20%. The effect is quite large, introduc- 
ing an error of ~ 0.15 in the determination of the reduced 
cosmological constant, and of ~ 1.5 in the determination of 
the global minimum for the equation of state wx (when no 
prior is assumed). 



7 DISCUSSION 

We discussed a method to probe the cosmological parame- 
ters, based on observations of giant luminous arcs in strong 
lensing clusters whose mass scales and mass profiles are de- 
termined with X-ray spectroscopic and surface brightness 
analyses. 

Giant arcs are more likely to form in massive, hot clus- 
ters with a small core radius. A large contribution to the 
energy budget of the universe in the form of dark energy 
with strong negative pressure also helps the occurrence of 
phenomena of strong lensing events. 

In order to constrain the cosmological parameters, we 
considered an updated subsample of rich, strong lensing 
clusters from the list in IWd i200Cf) . Only dynamically re- 
laxed, regular systems were included. Our statistical anal- 
ysis supports a dark energy dominated model of universe, 
in agreement with combined measurement from supernovae, 
galaxy clustering and CMB data, and a clear signature for an 
accelerated expansion of the universe is provided at the 2-(j 
confidence level. Dark energy with a strong negative pres- 
sure is also favoured. The cosmological constant provides the 
strongest candidate, but a wide range of quintessence mod- 
els are compatible. The data are consistent with the NEC 
at 3-(T statistical significance. 

The method that we have proposed is independent of 
the supernovae data and provides an alternative tool to 
probe the expansion of the universe. Our analysis is based 
on a very limited number of hypotheses and is not sensitive 
to Hq. Clusters of galaxies are not required to be standard 
rulers or candles. On the contrary of supernovae analyses, 
our analysis does not require either homogeneity in the sam- 
ple nor absence of evolution. Moreover, contamination from 
neighbours is truly negligible. 

Then, the redshift range covered by observations of 
lensed, highly magnified galaxies extends well beyond the 
supernova survey limits and can probe the changing point 
between the pressureless matter dominated era and the dark 
energy one. 

The method however needs a detailed modeling of the 
lens from X-ray observations and a deep understanding of 
the structure of the X-ray potential well is crucial. The iden- 
tification of irregularities and substructures in a complex 
morphology is very important when predicting arc position. 

Assumptions made in the methodology must be ver- 
ified testing systematics. The use of new X-ray data en- 
ables us to refine the criteria to select clusters in hydro- 
static equilibrium, such as the nearly coincidence of a X- 
ray single peak with the cD galaxy position or, mainly, the 
presence of a cooling-flow. Thus far, detailed X-ray analy- 
ses, based on Chandra or XMM-Newton observations, have 



been explicitly derived for several lensing clusters. Whereas 
for some clusters (A1835, MS 1358, and A239 0), there is 
a good agreement with the lensing analyses (jAllen et alJ 
[gpOla; Schinidt et al. 2001; Arabadiis et al..200a). a signifi- 
cant discre pancy is found, for exam ple, in the i rregular clus- 
ters A 2218 jMachacek et alj|2002h and A1689 JXue and Wul 
1200.^ . Whereas A2218 shows an offset between the X-ray 
centroid and the central dominant cD galaxy in A1689 the 
X-ray centre pe rfectly coincides with th e cD galaxy. An anal- 
ysis of AC 114 jPe Filippis et al]l2004) shows a good agree- 
ment between X-ray and lensing results, despite an offset 
of ~ 10". In a forthcoming paper, we want to extend our 
method to clusters of galaxies with an irregular X-ray mor- 
phology but very rich in multiple image systems. 
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